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1.0  SUMMARY 


1.1  Contract  Description 

This  contract  studies  unique  optical  computing  concepts  which  use  nonlinear 
optical  phenomena  to  perform  matrix  multiplication  and  to  provide  reconfigurable  opti¬ 
cal  interconnection.  The  study  focuses  on  the  use  of  real-time  holography  in  nonlinear 
media  such  as  photorefractive  crystals  for  optical  computing. 

1.2  Scientific  Problem 

By  incorporating  the  parallel  nature  of  optics  in  nonlinear  media,  it  is  possible 
to  perform  parallel  matrix  multiplication  using  four-wave  mixing.  In  addition,  the 
dynamic  holography  in  nonlinear  optical  media  provides  a  nature  candidate  for  the  recon¬ 
figurable  interconnection.  The  general  problem  in  this  program  is  to  generate  and  inves¬ 
tigate  new  concepts  which  use  these  nonlinear  optical  phenomena  for  optical  computing. 

Specifically,  this  program  investigates  experimentally  and  theoretically  the 
multiplication  of  matrices  using  optical  four-wave  mixing  in  nonlinear  media,  and  the 
possibility  of  using  such  matrix  multiplication  for  reconfigurable  interconnection. 

1.3  Progress 

There  are  several  areas  of  significant  progress  in  the  first  year  of  this  program 
that  are  directly  related  to  the  development  of  optical  matrix  multiplication  and  recon¬ 
figurable  interconnection.  These  include: 

1.  First  experimental  demonstration  of  parallel  matrix-vector  multiplier 
using  optical  four-wave  mixing  in  a  barium  titanate  (BaTiO^)  crystal. 

2.  First  experimental  demonstration  of  summation  process  inside  the  non¬ 
linear  media  in  matrix-vector  multiplications. 

3.  Experimental  demonstration  of  2  x  2  matrix-matrix  multiplication  using 
four-wave  mixing  in  a  BaTiO-j  crystal. 
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4.  Experimental  demonstration  of  matrix-matrix  multiplication  using  color 
multiplexing. 

5.  Experimental  demonstration  of  matrix-matrix  multiplication  using  convo¬ 
lution. 

6.  Development  of  a  new  scheme  of  reconfigurable  optical  interconnection 
which  is  energy-efficient. 

1.4  Publications  and  Presentations 


"Optical  Matrix-Vector  Multiplication  via  Four-Wave  Mixing  in  Photorefractive  Media," 
Opt.  Lett.  12,  138  (1987);  Opt.  Lett.  _12,  373  (1987). 

"Optical  Matrix-Vector  Multiplication  Using  Four-Wave  Mixing,"  Paper  MM5,  OSA 
Annual  Meeting  (Oct.  1986,  Seattle,  WA),  3.  Opt.  Soc.  Am.  A3  (13),  16  (1986). 

"Optical  Matrix-Matrix  Multiplication  Using  Multicolor  Four-W'ave  Mixing,"  Proc.  SPIE, 
vol.  881,  paper  38,  to  be  presented  at  OE/LASE’88  (Jan.  10-17,  1988,  Los  Angeles,  CA). 

1.5  Patent  Disclosures 

"Multicolor  Matrix-Matrix  Multiplier  with  Parallelism,"  M.  Khoshnevisan,  A.E.T. 
Chiou  and  P.  Yeh,  Rockwell  Patent  Disclosure  87SC37. 

"Matrix-Vector  Multiplier  using  Photorefractive  Phase  Conjugators,"  P.  Yeh,  A.E.T. 
Chiou  and  M.  Khoshnevisan,  Rockwell  Patent  Disclosure  87SC46. 

3 

"Optical  Matrix-Matrix  Multiplier  with  N  Parallelism  by  Spatial  Convolution  via  Four- 
Wave  Mixing,"  A.E.T.  Chiou,  Rockwell  Patent  Disclosure  87SC63. 

"Reconfigurable  Optical  Interconnect  using  Dynamic  Holograms,"  P.  Yeh,  Rockwell 
Patent  Disclosure  87SC54. 
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2.0  TECHNICAL  DISCUSSION 

The  technical  problem  addressed  in  tnis  contract  is  the  study  of  optical  matrix 
multiplication  and  reconfigurable  interconnection.  By  incorporating  the  parallel  nature 
of  optics  in  nonlinear  media  such  as  photorefractive  crystals,  ’  it  is  possible  to  perform 
parallel  matrix  multiplication  using  four-wave  mixing.  In  addition,  the  dynamic  holog¬ 
raphy  in  nonlinear  optical  media  provides  a  nature  candidate  for  the  reconfigurable  inter- 
3  4 

connection. 

Specifically,  this  program  is  aimed  at  the  utilization  of  nonlinear  optical 

5-& 

phenomena  such  as  phase  conjugation  and  four-wave  mixing  in  nonlinear  media  to 
perform  the  matrix-vector  and  matrix-matrix  multiplication.  Photorefractive  crystals 
such  as  BaTiO^  and  strontium  barium  niobate  (SBN)  were  used  as  the  nonlinear  media. 

Generally  speaking,  four-wave  mixing  is  a  nonlinear  optical  process  in  which 
three  input  waves  mix  to  yield  a  fourth  wave.  In  phase-matched  four-wave  mixing,  the 
three  input  waves  consist  of  two  counterpropagating  pump  waves,  Ej  and  E2,  and  an  arbi¬ 
trary  prohe  wave,  E3.  All  three  couple  through  the  third-order  susceptibility,  to 
yield  a  fourth  wave,  E^,  which  is  proportional  to  the  product  of  Ej,E2  and  complex  conju¬ 
gate  of  E3,  and  can  be  written  as^ 

E4  -  x(3)  E]E2E3* 

Such  a  four-w'ave  mixing  can  be  understood  in  terms  of  the  recording  and  readout  proces¬ 
ses  which  occur  in  holography.  In  the  nonlinear  media,  one  of  the  pump  waves  and  the 
probe  wave  form  an  interference  pattern  which,  in  turn,  induces  an  index  grating.  The 
other  pump  wave  is  diffracted  from  this  grating  and  generates  the  fourth  wave.  The 
formation  and  readout  processes  take  place  at  the  same  time.  Thus,  four-wave  mixing  is 
sometimes  referred  to  as  a  real-time  holography. 

By  using  four-wave  mixing  in  nonlinear  media,  multiplication  of  signals  can 
take  place  in  a  subpicosecond  time  scale.  In  addition,  if  we  use  the  parallel  nature  of 
optical  waves,  each  wave  can  carry  spatial  information  for  the  purpose  of  image  process¬ 
ing  using  the  two  transverse  dimensions.  Such  an  application  of  four-wave  mixing  to 
real-time  image  processing  was  recently  reported  using  BSO  crystals.  In  what  follows, 
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we  will  describe  some  unique  concepts  which  use  the  two  transverse  dimensions  to  carry 
the  matrix  information  for  the  purpose  of  matrix  multiplication  and  reconfigurable 
interconnection. 

2.1  Matrix-Matrix  Multiplication 

The  matrix  multiplication^"^  between  two  N  x  N  matrices  can  be  stated  as 

follows: 


C  =  AB 


where 


cij  ■  j  Bkj  <2> 

Note  that  a  matrix  multiplication  consists  of  two  main  operations,  a  parallel  multiplica¬ 
tion  and  an  addition. 

Referring  to  Fig.  la,  let  us  consider  a  four-wave  mixing  configuration  which  is 
suitable  for  matrix  multiplication.  Beam  1  and  Beam  2  contain  the  information  about  the 
two  matrices  A(x,z)  and  B(z,y),  respectively.  Beam  1  is  propagating  along  the  y-axis,  and 
Beam  2  is  propagating  along  the  x-axis.  These  directions  are  chosen  for  the  sake  of  clar¬ 
ity  in  introducing  the  concept.  They  are  not  the  only  directions  for  matrix  multiplica¬ 
tion.  Also,  these  matrices  can  be  either  continuous  or  discrete.  In  the  discrete  case, 
each  beam  consists  of  a  mattix  of  beamlets,  as  shown  in  Fig.  lb. 

In  the  nonlinear  medium,  these  two  matrix-carrying  beams  form  an  interfer¬ 
ence  pattern.  As  a  result  of  the  nonlinear  response  of  the  medium,  a  volume  grating  is 
formed.  This  grating  contains  information  about  the  product  of  the  matrix  elements  of 
these  two  matrices,  and  can  be  written  as: 


=  A(x,z)B*(z,y)ei +  c.c. 


C9028D/sn 
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where  Ris  the  difference  of  the  wavevector  of  the  matrix-carrying  beams  and  cc  repre¬ 
sents  complex  conjugation.  The  parameter  r\2  is  the  Kerr  coefficient  and  is  proportional 

(3) 

to  the  third-order  susceptibility  x'  of  the  medium.  Note  that  the  nonlinear  response  of 
the  medium  serves  the  purpose  of  parallel  multiplication. 

The  volume  grating  is  then  read  out  by  a  third  beam  which  can  simply  be  a 
plane  wave.  The  diffracted  beam  consists  of  the  integrated  contribution  from  each  part 
of  the  grating  along  the  beam  path,  and  thus  can  be  written 

C(x,y )  a  /  A(x,z)B*(z,y)dz  (4) 

where  the  integration  is  carried  out  along  the  beam  path.  Note  that  the  integration  com¬ 
pletes  the  operation  of  matrix  multiplication.  The  information  about  the  product  of 
these  two  matrices  is  now  impressed  on  the  transverse  spatial  distribution  of  the  dif¬ 
fracted  beam. 

Due  to  the  phase  matching  requirement,  the  readout  beam  must  be  incident 
along  the  directions  which  satisfy  the  Bragg  diffraction  ..undition  to  achieve  high  effi¬ 
ciency.  In  anisotropic  nonlinear  media,  the  polarization  states,  as  well  as  the  direction 
of  propagation,  can  be  chosen  such  that  the  largest  of  the  nonlinear  susceptibilities  is 
fully  utilized. 

The  four-wave  mixing  can  be  either  degenerate  or  nondegenerate.  In  the 
degenerate  case,  all  the  beams  have  the  same  frequency.  In  the  nondegenerate  case,  the 
frequencies  of  the  beams  can  be  slightly  different.  This  may  be  useful  for  the  purpose  of 
separating  the  diffracted  beam  from  the  undiffracted  portion. 

To  illustrate  the  information  capacity  of  such  a  nonlinear  optical  matrix 
multiplication,  let  us  consider  a  four-wave  mixing  process  using  an  Ar  ion  laser  at  4880 A 
in  a  medium  of  a  1  cm  cube.  The  grating  space  is  of  the  order  of  0.5  pm.  Thus,  10  pm  x 
10  pm  is  enough  for  each  pixel  of  information.  In  other  words,  a  1  cm  cube  is  cable  of 
handling  1000  x  1000  matrices.  With  a  material  response  time  of  1  ns,  such  a  matrix 
multiplication  has  a  potential  data  throughput  rate  of  quadrillion  bits  per  second 
(1015  bits/s)! 
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2.2  Matrix-Vector  Multiplication 

Another  scheme  as  shown  in  Fig.  2  is  suitable  for  matrix-vector  multiplica¬ 
tion.^”*^  Here,  as  an  example,  let  us  consider  a  discrete  case  in  which  we  need  to  carry 
out  the  multiplication  of  an  N-element  vector  and  an  N  x  N  matrix.  The  vector  is  fanned 
out  into  N-rows  of  identical  vectors.  These  NxN  small  beams  are  directed  to  a  non¬ 
linear  medium.  The  matrix  which  also  contains  NxN  small  beams  is  also  directed  to  the 
medium  in  such  a  way  that  each  beam  of  the  matrix  is  counterpropagating  in  a  direction 
relative  to  the  corresponding  beam  of  the  vector.  Thus,  in  the  medium,  there  are  N  x  N 
spatially  separated  regions  which  are  pumped  each  by  a  pair  of  counterpropagating 
beams.  Now,  N  x  N  probing  beams  are  directed  into  the  medium  in  such  a  way  that  each 
probe  beam  will  propagate  through  an  intersection  region.  The  probe  beams  will  be 
"plane  wave"  beamlets  propagating  in  parallel.  As  a  result  of  the  four-wave  mixing,  each 
probe  beam  will  generate  a  phase-conjugated  beam  which,  within  a  proportional  factor, 
can  be  written  M(i,j)a(j).  By  using  a  cylindrical  lens,  a  summation  over  j  can  be  ob¬ 
tained.  Thus,  we  have 

b(i)  =  l  M ( i , j ) a ( j )  (5) 

J 

where  a(j)  is  the  j-th  element  of  the  vector  a  and  M(i,j)  is  the  matrix  element.  Such  a 
scheme  for  matrix-vector  multiplication  can  also  be  used  for  matrix-matrix  multiplica¬ 
tion  by  decomposing  a  matrix  into  column  vectors  and  then  multiplying  the  matrix  with 
each  of  the  column  vectors. 

2.3  Recon figurable  Optical  Interconnection 


Optical  interconnection  will  play  a  key  role  in  both  the  optical  computing  and 
VLSI  systems. ^  There  are  many  advantages  of  optical  interconnect.  Optical 
systems  can  have  high  space-bandwidth  and  time-bandwidth  products;  hence,  many  inde¬ 
pendent  channels  that  could  be  exploited  for  demanding  computations.  Optical  proces¬ 
sors  are  inherently  two-dimensional  and  parallel.  Optical  signals  can  propagate  through 
each  other  in  separate  channels  with  essentially  no  interaction.  Optical  signals  can  also 
propagate  in  parallel  channels  without  any  interference  and  crosstalk.  In  the  VLSI 
systems,  optical  interconnect  can  be  used  to  solve  the  problem  of  communication,  as  well 
as  th‘  lock  distribution. 
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Fig.  2  Matrix-vector  multiplication  via  four-wave  mixing. 

Computer-generated  holograms  (CGH)  can  be  used  for  fixed  optical  intercon¬ 
nect.  The  most  general  CGH  consists  of  a  two-dimensional  array  of  subholograms.  Each 
subhologram  is  capable  of  diffracting  a  gate  output  to  any  gate  or  combination  of  gate 
inputs.  Such  a  fixed  interconnect  pattern  is  adequate  for  many  applications,  including 
the  point  operation,  matrix  operation,  the  clock  distribution  and  a  globally  synchronous 
systolic  processor.  However,  there  are  many  situations  in  which  a  dynamic  optical  inter¬ 
connect  is  required. 

In  Fourier  transform  and  sorting,  every  element  of  an  output  array  is  affected 
by  all  the  elements  of  the  input  array,  and  conversely,  each  element  of  the  input  array 
affects  all  elements  of  the  output  array.  In  addition,  the  interconnections  change  at  dif¬ 
ferent  stages  of  the  computation.  Thus,  these  operations  require  global  and  dynamic 
interconnection  between  different  elements  of  the  input  array. 
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In  image  restoration  and  pattern  recognition,  the  interconnect  patterns  could 
be  data-dependent,  making  it  impossible  to  foresee  the  interconnect  requirements  at 
different  stages  of  processing  without  having  foreknowledge  of  the  input.  The  computa¬ 
tional  throughput  of  parallel  processor  implementing  these  types  of  operation  will  be 
critically  affected  by  the  availability  of  a  dynamic  and  global  interconnect  network. 

In  Fourier  transform  and  sorting,  the  requirement  for  dynamic  interconnect 
can  be  avoided  by  resorting  to  a  fixed  but  global  interconnect  pattern  known  as  the  "per¬ 
fect  shuffle".  In  the  area  of  image  restoration  and  pattern  recognition,  there  is  a  mini¬ 
mum  amount  of  regularity  and  structure.  In  addition,  there  is  a  possible  data  and  time 
dependency  in  the  interconnect  requirements.  A  global  and  reconfigurable  interconnect 
network  will  be  vita!  to  achieving  high  throughput  and  high  efficiency  with  parallel  proc¬ 
essors  implementing  these  operations. 


The  most  general  interconnect  system  is  one  in  which  any  gate  output  can  be 
connected  to  the  input  of  any  gate  or  combination  of  gates  (see  Fig.  3).  The  effect  of 
such  an  interconnect  can  be  represented  by  the  matrix  equation 

O  =  M  I  (b) 

where  I  is  a  vector  representing  the  two-dimensional  input  array,  M  is  the  matrix  repre¬ 
senting  the  interconnect,  and  O  is  a  vector  representing  the  output  array.  In  digital  opt i - 
14  15 

cal  computing,  ’  the  input  array  is  actually  the  gate  output  array,  bach  matrix 
element  (i,j)  is  nonzero  if,  and  only  if,  there  is  a  connection  between  pixel  j  of  the  input 
array  and  pixel  i  of  the  output  array.  The  matrix-vector  multiplication  technique 
described  in  Section  2.2  can  be  used  as  an  optical  interconnect,  which  provides  both  local 
and  global  communication  between  the  gate  input  and  gate  output.  In  the  VLSI  systems, 
the  vectors  I  may  consist  of  N  laser  beams,  each  containing  a  stream  of  data  from  a 
processor,  and  the  vector  O  represents  the  output  laser  beams,  each  feeding  into  a  proc¬ 
essor.  The  main  advantage  of  using  such  nonlinear  matrix-vector  multiplication  for 
interconnection  is  that  the  matrix  M  can  be  easily  changed  for  the  purpose  of  reconfigur¬ 
ing  the  interconnect. 
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Fig.  3  Reconf  igurable  optical  interconnection. 


Equation  (6)  can  be  generalized  such  that  both  the  input  array  I  and  the  output 
array  O  are  two-dimensional  matrices.  In  that  case,  the  most  general  interconnect 
matrix  M  must  be  four-dimensional  (i.e.,  tensor  of  rank  four).  Each  matrix  element  (ijkl) 
represents  the  connection  between  pixel  (ij)  of  the  output  array  and  pixel  (kl)  of  the  input 
array.  Because  of  the  two-dimensional  nature  of  optical  waves,  it  is  desirable  to  repre¬ 
sent  the  matrix  M  by  a  two-dimensional  array  of  two-dimensional  submatrices.  Thus,  the 
matrix-matrix  multiplication  described  in  Section  2.1  can  be  used  for  reconfigurable 
optical  interconnect  which  provides  the  most  general  interconnection  pattern.  While  this 
interconnection  scheme  allows  complete  generality,  a  price  is  paid  in  terms  of  the  space- 
bandwidth  product  requirements  on  the  nonlinear  media. 
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3.0  PROGRESS 

3.1  Matrix-Vector  Multiplication 

In  our  earlier  demonstration  of  optical  matrix-vector  multiplication  (see 
Appendix  5.1),  we  have  used  optical  phase  conjugation  via  nonlinear  four-wave  mixing  in 
a  photorefractive  BaTiO^  crystal  to  perform  pixel-by-pixel  multiplication.  The  summa¬ 
tion  required  to  obtain  matrix-vector  products  is  performed  subsequently  by  a  cylindrical 
lens  external  to  the  nonlinear  medium.  With  a  slight  modification  of  the  experimental 
geometry,  we  have  successfully  demonstrated  that  both  the  pixel-by-pixel  product  and 
the  summation  can  be  done  inside  the  nonlinear  medium  to  achieve  the  desired  matrix- 
vector  product  without  an  external  cylindrical  lens. 

The  top  view  of  the  new  experimental  geometry  is  illustrated  schematically  in 
Fig.  4.  The  readout  vector  in  this  case  consists  of  a  vertical  column  of  beamlets  with 
equal  intensity  and  enter  the  crystal  at  an  oblique  angle  such  that  each  beamlet 
traverses,  inside  the  crystal,  all  the  counterpropagating  pumping  beamlets  at  the  same 


PHOTOREFRACTIVE 


Fig.  4  Optical  matrix-vector  multiplication  using  photorefractive  four-wave  mixing. 
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horizontal  plane.  As  the  proper  elements  of  the  matrix  and  the  vector  are  encoded  in  the 
two  counterpropagating  beams  via  the  appropriate  masks,  the  phase-conjugate  output  of 
each  beamlet  consists  of  the  sum  of  the  product  resulting  from  each  individual  encoun¬ 
ter.  From  a  pure  geometrical  point  of  view,  the  probe  beam  (i.e.,  the  readout  vector) 
can  conveniently  be  injected  (from  c-face  of  the  crystal)  perpendicular  to  the  counter- 
propagating  pumping  beams.  The  tensor  nature  of  the  photorefractive  grating, however, 
does  not  allow  coupling  in  such  a  geometry,  except  for  cubic  crystals  such  as  GaAs. 

The  experimental  result  for  the  multiplication  of  a  A  x  4  binary  matrix  and  a 

4  x  1  binary  vector  is  shown  in  Fig.  5.  The  magnified  images  of  the  masks  used  to  encode 

the  matrix,  the  vector  and  the  probe  are  shown  in  Figs.  5a,  b  and  c,  respectively.  The 

experimental  results  representing  the  final  matrix-vector  product  is  shown  in  Fig.  5d. 

The  relative  intensity  of  each  output  spot  representing  the  element  of  the  product  vector 

is  shown  in  Fig.  6.  Due  to  the  slow  response  time  (or  the  order  of  second  for  optical 

2 

intensity  of  the  order  of  a  few  tens  of  milliwatt/cm  )  and  the  consequential  sensitivity  to 
environmental  changes  of  the  photorefractive  process,  the  output  intensity  fluctuates 
significantly.  The  ±25%  noise  level  is  likely  to  limit  its  application  to  binary  operation 
only. 

3.2  Matrix -Matrix  Multiplication  via  Color  Multiplexing 

We  have  also  extended  our  experimental  demonstration  on  matrix-vector 
multiplication  to  matrix-matrix  multiplication  by  color  multiplexing  (see  Appendix  5.2). 
The  basic  idea  is  to  decompose  the  problem  into  matrix-vector  multiplications,  as  shown 
in  Fig.  7,  and  carry  out  all  the  matrix-vector  multiplications  simultaneously  in  parallel  by 
using  color  multiplexing.  Note  that  each  of  the  matrix-matrix  operations  shown  on  the 
right-hand  side  of  the  equation  in  Fig.  7  is  in  fact  a  matrix-vector  multiplication.  The 
basic  principle  of  color  multiplexing  used  to  encode  the  component  vectors  with  different 
colors  is  illustrated  in  Fig.  8  for  the  case  of  4  x  4  matrices.  Each  row  vector  of  the 
matrix  Mj  is  illuminated  by  one  color,  and  all  the  color  components  are  then  combined  by 
the  prism  (angular  multiplexing)  into  a  single  row  prior  to  further  expansion  by  ana- 
morphic  optics  (not  shown  in  the  figure)  to  match  the  mask  representing  the  second 
matrix.  After  proper  element-by-element  multiplications  and  summation  (summing 
optics  omitted)  into  the  column,  the  resulting  multicolor  output  is  demultiplexed  into  dif¬ 
ferent  color  components  that  together  represent  the  final  product. 
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M3  =  MtM2 


Fig.  8  Optical  matrix-n.c:  ix  multiplication  via  color  multiplexing. 

Using  an  argon  ion  laser  that  oscillates  simultaneously  at  five  colors  in  the 
blue-green  and  a  BaTiO-j  crystal,  we  have  demonstrated  the  principle  described  above  for 
the  case  of  2  x  2  matrices.  The  two  matrices  M  and  A  and  the  corresponding  color 
encoding  schemes  are  shown  in  Fig.  9,  and  the  experimental  results  are  given  in  Fig.  10. 
The  spatial  pattern  in  Fig.  10a  is  the  output  prior  to  the  summing  cylindrical  optics;  the 
pattern  in  Fig.  10b  is  the  final  result.  The  relative  intensity  of  each  spot  is  also  shown  in 
the  figure. 

An  important  practical  issue  is  the  maximum  dimension  (N)  of  the  matrices 
allowed  for  this  approach.  For  N  <  5,  photorefractive  materials  such  as  BaTiC>3  or  SBN 
can  be  used.  For  larger  N,  however,  it  is  expected  that  the  diffraction  efficiency  of  the 
hologram  for  any  given  color  will  be  adversely  affected  by  the  erasure  effects  of  the 
other  N-l  beams  and  by  the  superposition  of  the  other  N-l  holograms.  For  larger  values 
of  N,  modification  such  as  the  use  of  near-resonance  multilayer  nonlinear  medium  (mul¬ 
tiple  quantum  well,  for  instance)  may  be  helpful. 
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Fig.  9  The  color  encoding  scheme  used  in  our  experimental  demonstration  of  optical 
matrix-matrix  multiplication  for  the  case  of  2  x  2  matrices. 
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Lxperimental  results  for  matrix-matrix  multiplication  using  color-multiplexed 
four-wave  mixing  in  a  photorefractive  crystal:  (a)  the  spatial  intensity  pattern 
of  the  output  prior  to  the  scanning  cylindrical  optics;  (b)  the  spatial  pattern 
representing  the  final  result. 
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3.3  Matrix -Matrix  Multiplication  using  Spatial  Convolution 

In  addition  to  the  matrix-matrix  multiplication  using  color-multiplexed  four- 
wave  mixing  described  in  the  previous  section,  we  have  also  successfully  demonstrated 
another  scheme  for  matrix-matrix  multiplication  using  spatial  convolution  via  four-wave 
mixing.  From  the  experimental  point  of  view,  the  key  difference  between  this  approach 
and  the  others  described  above  is  that  the  nonlinear  crystal  is  now  located  at  the  com¬ 
mon  Fourier  plane,  rather  than  the  common  image  plane,  of  the  input  matrix  masks.  The 
encoding  scheme,  as  explained  in  the  following  paragraphs,  is  also  different  from  those 
used  in  the  earlier  approaches.  From  the  conceptual  point  of  view,  matrix-matrix  multi¬ 
plication  in  full  parallelism  is  achieved  by  space-multiplexing  via  spatial  convolution 
using  degenerate  four-wave  mixing. 

White  and  Yariv^  have  demonstrated  that  (spatial)  convolution  and 
correlation  of  two  (two-dimensional)  patterns  can  be  achieved  in  real  time  by  four-wave 
mixing  in  the  common  Fourier  plane  of  the  input  patterns  and  recording  the  phase- 
conjugated  output  at  the  corresponding  object  plane.  A  typical  experimental 
configuration  is  shown  schematically  in  Fig.  11.  c  '-'cifically,  if  U3  is  a  small  aperture 
simulating  a  point  source  (or  a  delta  function),  a  pattern  representing  the  convolution  of 
U1  and  U2  (U1*U2)  is  observed  in  the  outptt  plane.  An  illustrative  example  is  given  in 
Fig.  12.  The  design  of  U1  and  U2  to  represent  two  matrices  so  that  their  product  is 
represented  by  U1*U2  is  explained  in  the  following  paragraphs. 

For  simplicity  and  clarity,  let  us  consider  the  simplest  caseb  of  two  matrices  A 
and  B,  both  of  dimension  2x2,  and  their  product  C  =  AB  given  below. 


A  =  a. 


B  =  b, 


C  =  AB 


:11  c  1 2 


aj  ,b| ]  + a,2b2| 
a  2 1 b 1 1  +  a22b21 


a  1 1 b 1 2  +  a12b22 
a  2 1 b2 1  +  a22b22 
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-NONLINEAR  SC41299 

OPTICAL  MEDIUM 
(AT  THE  COMMON  FOURIER  PLANE) 


OUTPUT  PLANE 
(FOCAL  PLANE  OF  L3) 


A  schematic  diagram  illustrating  the  basic  idea  of  four-wave  mixing  in  a  non¬ 
linear  medium  located  at  the  common  Fourier  plane  of  the  input  spatial 
pattern. 


U1  *  U2 


A  schematic  illustration  of  convolution  of  two  spatial  patterns:  (a)  and  (b)  are 
the  two  input  spatial  patterns;  (c)  represents  the  resulting  spatial  convolution 
of  the  two  inputs. 
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A  transparency  Li  1  consists  of  four  small  apertures,  each  with  its  intensity  transmittance 
proportional  to  each  of  the  matrix  element  a  is  shown  in  Fig.  13a.  A  similar  trans¬ 
parency  L'2  corresponding  to  the  transpose  of  B  (i.e.,  rows  and  columns  interchanged)  is 
shown  in  Fig.  13b.  Note  that  the  vertical  distance  between  the  elements  in  U1  is 
designed  so  that  it  is  considerably  larger  than  that  in  U2,  while  the  horizontal  distances 
between  the  elements  are  identical  in  U1  and  U2.  Note  also  that  the  dashed  rectangular 
box  connecting  the  four  elements  in  U2  is  an  artifact  to  facilitate  the  following 
explanation  and  should  be  ignored  in  the  actual  design.  By  comparing  Figs.  12  and  13,  it 
is  straightforward  to  see  how  the  quadruplet  of  the  dashed  rectangular  box  is  formed  at 
the  corresponding  position  in  Fig.  13c  as  a  result  of  the  two-dimensional  spatial 
convolution.  As  noted  before,  all  the  dashed  lines  in  Fig.  13c  are  artifacts  and  the  actual 
result  should  show  only  the  corners  of  each  box,  i.e.,  a  total  of  16  spots  with  the  intensity 
of  each  proportional  to  the  cross  product  of  all  the  elements  in  A  and  those  in  B.  The 
four  doublets  lying  along  the  Y-axis  in  Fig.  13c  are  drawn  slightly  offset  from  the  axis  to 
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Fig. 


13  A  schematic  diagram  illustrating  optical  matrix-matrix  multiplication  bv 
convolution  for  the  case  of  2  x  2  matrices:  (a)  matrix  A  and  the  mask  U  ,  to 
encode  A;  (b)  matrix  B  and  the  mask  U?  to  encode  B  ;  (c)  the  spatial  pattern 
resulting  from  convolution  of  Uj  and  U?. 
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expose  the  individual  component.  In  practice,  the  two  components  of  each  doublet  are 
spatially  overlapping  on  the  Y-axis.  The  intensity  of  the  four  doublets,  from  top  to 
bottom,  on  the  Y-axis  are  proportional  to  c12,  Cjj,  c22  and  c21,  respectively.  The  eight 
cross  terms  that  do  not  contribute  to  the  matrix-matrix  multiplication  can  be  filtered 
out  easily  as  they  are  physically  separated  from  the  Y-axis. 

Another  example  illustrating  the  case  of  multiplying  a  2  x  3  matrix  with  a 
3x4  matrix  is  shown  in  Fig.  14.  Note  that  all  the  cross  terms  that  do  not  contribute  to 
the  matrix-matrix  product  are  omitted  in  Fig.  14c  for  the  sake  of  clarity. 

Using  an  argon  ion  laser  (5145A)  and  a  BaTiC>3  crystal,  we  have  experimentally 
demonstrated  the  concept  described  above  for  the  case  of  2  x  2  matrices.  The  images  of 
the  masks  used  to  encode  the  two  matrices  A  and  B  are  shown  in  Figs.  15a  and  15b, 
respectively.  The  experimental  result  representing  the  product  is  shown  in  Fig.  i5c.  An 
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Fig.  14  An  illustration  similar  to  Fig.  13  for  the  case  when  the  dimensions  of  the 
matrices  A,  B  and  C  are  2  x  3,  3  x  4  and  2x4,  respectively. 
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Fig.  15  Experimental  results  for  optical  matrix-matrix  multiplication  by  convolution 
via  four-wave  mixing  in  the  spatial  frequency  domain,  (a)  and  (b)  are  the 
matrix  masks  used  for  encoding  the  matrices  A  and  BT;  (c)  shows  the  result  of 
their  spatial  convolution;  (d)  is  an  overexposed  version  of  (c)  to  reveal  the 
cross  terms  that  do  not  contribute  to  the  matrix  product. 

overexposed  version  of  Fig.  15c  is  shown  in  Fig.  15d  to  reveal  the  noise  resulting  from 
the  other  cross  terms  that  do  not  contribute  to  the  matrix  product. 

In  principle,  the  basic  concept  described  above  can  be  appled  to  the  multiplica¬ 
tion  of  a  matrix  (of  any  dimension  mxn)  and  a  compatible  matrix  B  (of  dimension  nxp).  In 
practice,  the  signal-to-noise  (S/N)  ratio  is  expected  to  degrade  as  the  dimension  of  the 
matrix  increases.  The  problem  of  S/N  is  less  serious  in  the  mixed  binary  mode.  The 
major  advantage  of  this  scheme  is  the  full  parallelism  achieved  with  relativey  simple 
optical  arrangements.  From  the  technological  point  of  view,  the  fact  that  all  the  output 
elements  lie  on  a  straight  line  can  be  significant,  since  this  allows  one  to  use  a  linear 
detector  array  instead  of  a  2-D  array  for  output  detection. 


21 

C9028D/sn 


Rockwell  International 

Science  Center 


SC5502.AR 

3.4  Optical  Interconnection  using  Dynamic  Photorefractive  Holograms 

During  the  first  year  of  the  program,  we  have  also  developed  a  new  concept  of 
reconfigurable  optical  interconnection  using  photorefractive  holograms.  This  new 
scheme  provides  an  optical  interconnection  between  an  array  of  laser  sources  and  an 
array  of  detectors  with  a  very  high  energy  efficiency. 

Reconfigurable  interconnection  linking  laser  arrays  and  detector  arrays  plays  a 
key  role  in  optical  computing.  Conceptually,  such  interconnection  can  be  achieved  by 
using  an  optical  matrix-vector  multiplication. 

v'  =  M  v  (8) 

where  v  is  the  input  vector  representing  the  signals  carried  by  the  array  of  lasers,  and  v' 
is  the  output  vector  representing  the  signals  carried  by  the  array  of  detectors.  The 
matrix  M  represents  the  interconnection  pattern.  When  a  transparency  or  spatial  light 
modulator  (SLM)  is  used  as  the  interconnection  pattern,  a  large  fraction  of  energy  is 
absorbed  by  the  transparency  or  SLM.  This  energy  loss  increases  as  the  dimension  of  the 
array  increases.  Often,  the  fractional  energy  loss  can  be  as  large  as  (N-D/N,  where  N  is 
the  dimension  of  the  array.  For  a  1000  x  1000  crossbar  switch,  the  loss  due  to  fanout  can 
be  as  big  as  99.9%.  This  is  not  acceptable  in  high-speed  computing  because  signals  are 
passing  through  the  SLM  billion  times  per  second  and  the  energy  loss  can  be  enormous. 

Referring  to  Figs.  16  and  17,  we  describe  a  new  method  of  reconfigurable  opti¬ 
cal  interconnection  which  uses  the  nonreciprocal  energy  transfer  in  photorefractive  two- 
wave  mixing  to  improve  the  energy  efficiency.  Figure  16  describes  a  one-dimensional 
case  for  the  sake  of  clarity  in  explaining  the  concept.  A  small  fraction  of  a  laser  beam  is 
coupled  out  of  the  beam  by  using  a  beam  splitter.  This  small  fraction  (called  probe 
beam)  is  then  expanded  by  using  a  cylindrical  lens  and  passes  through  the  spatial  SLM.  In 
the  example  shown,  the  laser  beam  is  to  connect  to  Detectors  b  and  d  as  prescribed  b> 
the  SLM.  The  transmitted  beam  is  then  recombined  with  the  main  beam  inside  a  photo¬ 
refractive  crystal.  As  a  result  of  nonreciprocal  energy  coupling,  almost  all  the  energy  in 
the  main  beam  is  transferred  to  the  probe  beam  which  carries  the  interconnection  pat¬ 
tern.  The  results  is  an  optical  interconnection  with  high  energy  efficiency.  Figure  17 
describes  the  reconfigurable  interconnection  for  laser  arrays  and  detector  arrays.  In  the 
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Fig.  16  A  schematic  drawing  of  a  1  x  N  optical  interconnection  using  dynamic 
photorefractive  holograms.  N  =  U. 
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Fig.  17  A  schematic  drawing  of  a  N  x  N  optical  interconnection  using  dynamic 
photorefractive  holograms.  N  =  4. 
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example  (a  4  x  4  interconnection)  shown,  Laser  1  is  to  connect  to  Detectors  b  and  c, 
Laser  2  is  to  connect  to  Detectors  a  and  d,  Laser  3  is  to  connect  to  Detectors  c  and  d, 
Laser  4  is  to  connect  to  Detector  a  and  c.  A  cylindrical  lens  is  used  to  focus  the  two- 
dimensional  array  of  beams  into  a  vector  (one-dimensional  array).  As  a  result,  Detector 
a  receives  signals  front  Lasers  2  and  4,  Detector  b  receives  signals  from  Laser  1,  Detec¬ 
tor  c  receives  signals  from  Lasers  1,  3  and  4,  and  Detector  d  receives  signals  from  Lasers 
2  and  3.  Such  a  concept  can  be  extended  to  interconnect  NL  lasers  with  ND  detectors 
where  NL  and  ND  are  two  large  numbers. 

The  interconnection  can  be  reconfigured  by  using  a  different  SLM  pattern. 
The  energy  loss  due  to  SI  is  no  more  than  the  reflectivity  of  the  beam  splitter  plus  the 
material  bulk  absorption  due  to  the  photorefractive  crystals.  The  beam  splitter  can  be 
chosen  such  that  the  reflectivity  is  small  (e.g.,  5%  or  less)  so  that  such  energy  loss  is 
minimized.  The  hologram  formation  time  will  limit  the  reconfiguration  time.  Once  the 
interconnection  pattern  is  formed  inside  the  photorefractive  crystal  as  a  hologram,  such 
a  scheme  is  capable  of  providing  the  interconnection  for  high  data  rate  transmission.  In 
such  an  interconnection,  the  output  of  each  laser  can  be  input  to  any  one  or  all  of  the 
detectors. 

Optical  phase  conjugation  can  also  be  used  in  conjunction  with  the  two-wave 
mixing  to  correct  for  any  phase  aberration  that  may  be  caused  by  the  crystal  imperfec- 
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in  Photorefractive  Media 
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We  propose  and  describe  a  new  method  of  optical  matrix-vector  multiplication  by  using  four-wave  mixing  in 
photorefractive  media.  Using  a  BaTiO.i  crystal,  we  have  demonstrated  such  a  parallel  multiplication.  The  results 
are  presented  and  discussed. 


A  large  number  of  signal-  and  image-processing  algo- 
rithi  s  can  be  expressed  in  terms  of  matrix  operations. 
The  multiplication  of  two  matrices  is  one  of  the  most 
basic  operations  in  matrix  algebra.  Optics,  v  h  its 
inherent  parallelism,  can  offer  great  improvement  in 
the  speed  of  these  operations.  Optical  processors  for 
multiplying  two  matrices  have  been  described  in  the 
literature.1-3  Considerable  work  has  also  been  report¬ 
ed  on  performing  optical  matrix-vector  multiplica¬ 
tion.3-5  In  this  Letter,  we  describe  a  new  method  of 
performing  matrix-vector  multiplication  by  using  op¬ 
tical  four-wave  mixing  in  nonlinear  media.  In  addi¬ 
tion,  we  also  report  the  first  experimental  demonstra¬ 
tion  to  our  knowledge  of  such  an  optical  matrix  proces¬ 
sor  by  using  a  photorefractive  BaTiCL  crystal. 

Optical  four-wave  mixing  has  been  a  subject  of  con¬ 
siderable  interest  during  the  past  several  years.  Much 
attention  has  been  focused  on  wave-front-correction6  ' 
and  phase-conjugate  interferometry.910  In  the  area  of 
signal  processing,  optical  four-wave  mixing  has  been 
used  to  perform  spatial  information  processing,11  im¬ 
age  subtraction,12-15  and  logic  operations.16  However, 
little  attention  has  been  paid  to  the  use  of  the  inher¬ 
ently  parallel  multiplication  nature  of  four-wave  mix¬ 
ing  for  optical  matrix  multiplication.  In  this  Letter 
we  describe  a  new  method  of  performing  matrix-vec¬ 
tor  multiplication  using  four-wave  mixing  in  nonlinear 
media. 

In  what  follows,  we  briefly  describe  the  basic  princi¬ 
ples  of  real-time  matrix  multiplication  using  optical 
four-wave  mixing  in  nonlinear  media.  For  the  sake  of 
convenience,  we  will  limit  ourselves  to  the  case  of 
square  matrices,  although  extension  to  nonsquare  ma¬ 
trices  will  be  straightforward. 

Referring  to  Fig.  1,  we  consider  a  scheme  that  is 
suitable  for  a  matrix-vector  multiplication.  Here,  as 
an  example,  let  us  consider  a  discrete  case  in  which  we 
need  to  carry  out  the  multiplication  of  an  N-element 
vector  and  an  N  X  N  matrix.  In  this  scheme,  the 
vector  is  fanned  out  into  N  rows  of  identical  vectors. 
These  N  X  N  beamlets  are  directed  to  a  nonlinear 
medium.  The  matrix,  which  also  contains  N  X  N 
beamlets,  is  also  directed  to  the  medium  in  such  a  way 
that  each  beamlet  of  the  matrix  is  counterpropagating 
in  a  direction  relative  to  the  corresponding  beamlet  of 


the  vectc'.  Thus,  in  the  medium,  there  are  N  X  N 
spatially  separated  regions,  each  of  which  is  pumped 
by  a  pair  of  counterpropagating  beamlets.  Now  a 
probing  beam,  which  consists  of  N  X  N  beamlets,  is 
directed  into  the  medium  in  such  a  way  that  each 
probe  beamlet  will  propagate  th.  ough  the  correspond¬ 
ing  intersection  region,  as  shown  in  Fig.  1.  The  probe 
beam  will  be  plane-wave  beamlets  propagating  in  par¬ 
allel.  All  the  beamlets  in  this  probe  beam  are  of  equal 
intensity.  As  a  result  of  the  four-wave  mixing,  each 
probe  beamlet  will  generate  a  phase-conjugated  beam- 
let,  which,  within  a  proportional  factor,  can  be  written 
as  M(i,  j)a(J).  By  using  a  cylindrical  lens,  a  summa¬ 
tion  over  j  can  be  obtained.  Thus  we  have 

6(0  =  V  M(i,j)  X  a(j),  (1) 

i 

where  a(j)  is  the  ;th  element  of  the  vector  a  and  M(i,  j) 
is  the  matrix  element.  Such  a  scheme  for  matrix- 
vector  multiplication  can  also  be  used  for  matrix-ma¬ 
trix  multiplication  by  decomposing  a  matrix  into  col¬ 
umn  vectors  and  then  multiplying  the  matrix  with 
each  of  the  column  vectors. 

The  probe  beam  can  also  simply  be  a  unif  -m  plane 
wave,  without  beamlets.  The  phase-conjugated  beam 
will  consist  of  N  X  N  beamlets,  because  only  a  matrix 
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Fig.  2.  Experimental  set  up  used  to  demonstrate  the  optical 
matrix-vector  multiplication. 


of  N  X  IV  regions  is  pumped  by  counterpropagating 
beamlets.  The  matrix  and  the  vector  information  can 
actually  be  carried  on  any  two  of  the  three  incident 
beams  in  the  four-wave  mixing  process. 

The  summation  process  can  also  be  obtained  with¬ 
out  the  external  cylindrical  lens  by  using  a  different 
scheme,  described  as  follows.  Consider  a  probe  beam 
that  consists  of  a  column  of  N  equal  beamlets.  The 
probe  is  incident  into  the  nonlinear  medium  in  such  a 
way  that  each  of  the  probe  beamlets  is  made  to  propa¬ 
gate  through  a  row  of  N  intersection  regions.  The 
phase-conjugation  process  automatically  performs  the 
summation  as  well  as  the  multiplication.  The  phase- 
conjugated  beam  is  thus  the  product  of  the  matrix- 
vector  multiplication. 

As  a  result  of  the  nature  of  the  four-wave  mixing 
process  in  the  medium,  this  new  matrix-vector  multi¬ 
plier  operates  on  the  field  amplitudes  and  thus  can  be 
used  to  handle  matrices  and  vectors  with  complex 
elements.  It  is  a  coherent  device  rather  than  an  inco¬ 
herent  one.  This  aspect  is  distinctly  different  from 
most  of  the  earlier  approaches,  which  are  all  incoher¬ 
ent.  When  the  device  is  operated  in  the  coherent 
mode,  the  phase  of  each  beamlet  must  be  maintained 
uniformly  over  the  transverse  dimension  of  the  beam- 
let.  In  addition,  such  phases  must  also  be  maintained 
fixed  in  the  summation  process.  If  these  phases  are 
not  uniform  over  the  beamlets,  the  final  step  becomes 
an  incoherent  summation  as  a  result  of  the  spatial 
averaging.  Under  such  circumstances,  this  matrix- 
vector  multiplier  operates  on  the  intensities  and  thus 
handles  only  positive  numbers. 

The  matrix-vector  multiplication  is  demonstrated 
experimentally  by  using  four-wave  mixing  in  a  photo- 
refractive  BaTiC>3  crystal.17  The  crystal  is  cut  in  such 
a  way  that  the  faces  are  all  perpendicular  to  the  princi¬ 
pal  axes.  All  the  beams  are  polarized  extraordinary  in 
the  xz  plane  and  are  incident  onto  the  a  faces  of  the 
crystal. 

In  our  experiments,  an  Ar-ion  laser  beam  at  wave¬ 
length  514.5  nm  with  an  output  power  of  a  few  hun¬ 
dred  milliwatts  is  used  as  the  light  source.  The  laser 
beam  is  expanded  and  collimated  into  a  beam  size  of 
approximately  1  cm.  Figure  2  shows  our  experimental 
setup.  The  expanded  beam  is  then  split  into  three 
beams  by  using  beam  splitters  BSi  and  BS2.  To  dem¬ 


onstrate  the  principle  of  operation,  we  chose  the  fol¬ 
lowing  matrix  and  vector: 
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The  matrix  and  the  vector  information  is  imprinted 
onto  two  of  the  beams  by  using  transparencies  (MA* 
and  M  A2  in  Fig-  2)  that  consist  of  circular  dots.  White 
circular  dots  represent  the  l’s,  and  the  dark  regions 
the  0’s.  Instead  of  using  a  cylindrical  lens  to  fan  out 
the  vector  into  five  identical  rows  as  described  earlier, 
we  simply  use  a  transparency  of  e  matrix  that  consists 
of  five  identical  rows  of  (1  0  1  1  1).  The  incident 
matrices,  corresponding  to  the  matrix  Af  and  the  vec¬ 
tor  a,  are  shown  in  Figs.  3(a)  and  3(b).  Four  lenses,  Li, 
L_>,  L3,  and  L4  (with  focal  lengths  f\  =  f3  -  0.5  m  and  /2 
=  f\  =  0.2  m)  are  used  to  image  the  transparencies  MAi 
and  MA,  onto  the  center  of  the  nonlinear  crystal, 
where  the  diameter  of  each  beamlet  is  0.18  mm  with  a 
center-to-center  distance  of  0.51  mm.  The  size  of  the 
whole  image  inside  the  crystal  is  2.2  mm  X  2.2  mm. 
The  uniform  beam  is  also  reduced  (and  collimated)  to 
a  diameter  of  3.6  mm  inside  the  crystal  by  lenses  L5 
and  L6  (with  focal  lengths  /5  =  0.6  m  and  fe  =  0.25  m, 
respectively).  The  angle  8  between  the  uniform  beam 
and  one  of  the  pump  beams  is  11°.  The  two  beams 
that  carry  the  matrix  and  the  vector  information,  re¬ 
spectively,  are  incident  onto  the  opposite  a  faces  of  a 
photorefractive  BaTiOa  crystal.  The  incident  angle  <t> 
is  approximately  20°  with  respect  to  the  surface  nor¬ 
mal.  As  a  result  of  the  four-wave  mixing,  a  phase- 
conjugated  beam,  consisting  of  a  5  X  5  beamlet  pat¬ 
tern,  is  generated  and  is  shown  in  Fig.  3(c).  A  cylin¬ 
drical  lens  is  then  used  to  pe  form  the  summation. 
The  resulting  product  vector  is  shown  in  Fig.  3(d).  By 


Fig.  3.  Photographs  showing  the  matrix  and  the  vectors: 
(a)  input  matrix;  (b)  input  vector,  expanded  into  a  matrix 
with  five  identical  rows;  (c)  resultant  phase-conjugate  reflec¬ 
tion  matrix;  (d)  resultant  vector  of  a  matri:  -vector  multipli¬ 
cation  after  summation  of  (c)  by  a  cylindrical  lens. 
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measuring  the  intensity  of  the  individual  dots,  we 
found  that  the  product  vector  is  (3  2  0  1  4),  in  good 
agreement  with  our  prediction;  in  other  words, 
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The  intensities  of  the  product  vector  elements  are 
measured  by  integrating  over  the  cross  section  of  each 
beamlet.  The  experimental  error  is  less  than  5%. 
This  small  error  is  due  to  environmental  perturbations 
and  the  slow  response  of  the  photorefractive  crystal. 
Although  our  experiment  demonstrates  the  operation 
for  matrices  and  vectors  with  binary  numbers  (1,  0), 
this  matrix-vector  multiplier  can  also  be  operated  in 
the  analog  mode  to  handle  elements  with  gray  levels. 
Such  an  analog  operation  in  photorefractive  crystals 
requires  that  the  probe-beam  intensity  be  weak  so  that 
the  error  due  to  the  intensity  denominator  in  the  index 
modulation  can  be  neglected.  This  material-related 
error  disappears  when  the  matrix  and  the  vector  con¬ 
tain  only  binary  numbers  (1,0). 

In  conclusion,  we  have  proposed  and  demonstrated 
a  new  method  of  matrix  operation  by  using  optical 
lour-wave  mixing  in  nonlinear  media.  In  particular, 
we  have  demonstrated  experimentally  matrix-vector 
multiplication  in  a  photorefractive  BaTiO,t  crystal. 

The  authors  acknowledge  helpful  discussions  with 
M.  D.  Ewbank  and  Ni.  khoshnevisan  (Rockwell  Inter¬ 
national)  and  J.  Neff  (Defense  Advanced  Research 
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SUMMARY 


A  new  method  for  multiplying  two  NxN  matrices  optically  with  N3  parallelism  is  described  and 


demonstrated  for  the  case  of  2x2  matrices  in  mixed  binary  modes.  The  basic  idea  is  to  decompose 


the  original  problem  into  N  problems  of  matrix-vector  multiplications  which  are  carried  out 


simultaneously  (in  parallel)  using  multi-color  four-wave  mixing.  Recently,  we  have  reported  a 


new  method  of  performing  optical  matrix-vector  multiplication  [1]  using  four-wave  mixing  in  a 


nonlinear  optical  medium.  For  media  which  have  significant  nonlinear  optical  susceptibility  over  a 
range  of  wavelengths  or  at  more  than  one  wavelength  (color),  multi-color  four-wave  mixing  [2] 
can  be  achieved;  and  several  matrix-vector  multiplications  can  be  carried  out  simultaneously.  We 
introduce  a  generic  scheme  wnich  can  extend  such  matrix-vector  multiplications  to  matrix-matrix 
multiplication.  Using  a  combination  of  angular  and  color  multiplexing/demultiplexing  techniques 
together  with  a  multi-color  four-wave  mixing  process,  we  show'  the  multiplication  of  two  matrices 
in  parallel  w'ith  time  complexity  of  the  order  of  N3.  Experimental  results  for  the  case  of  2x2 
matrices  using  tw-o-color  four-wave  mixing  in  a  photorefractive  barium  titanate  crystal  are  given. 


Some  of  the  difficulties  and  the  practical  limit  on  the  dimension  (N)  of  the  matrix  are  discussed. 
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